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ABSTRACT 

Neutron stars (NSs) in the astrophysical universe are often surrounded by accretion disks. Accretion 
of matter onto an NS may increase its mass above the maximum value allowed by its equation of 
state, inducing its collapse to a black hole (BH). Here we study this process for the first time, in 
three-dimensions, and in full general relativity By considering three initial NS configurations, each 
with and without a surrounding disk (of mass ~ 7% Mns), we investigate the effect of the accretion 
disk on the dynamics of the collapse and its imprint on both the gravitational wave (GW) and 
electromagnetic (EM) signals that can be emitted by these sources. We show in particular that, even 
if the GW signal is similar for the accretion induced collapse (AIC) and the collapse of an NS in 
vacuum (and detectable only for Galactic sources) , the EM counterpart could allow us to discriminate 
between these two types of events. In fact, our simulations show that, while the collapse of an NS 
in vacuum leaves no appreciable baryonic matter outside the event horizon, an AIC is followed by a 
phase of rapid accretion of the surviving disk onto the newly formed BH. The post-collapse accretion 
rates, on the order of ~ 1CU 2 Mq s _1 , make these events tantalizing candidates as engines of short 
gamma-ray bursts. 

Keywords: accretion, accretion disks — gamma-ray burst: general — gravitational waves — methods: 
numerical — stars: neutron 



1. INTRODUCTION 

In the astrophysical universe, neutron stars (NSs) are 
often surrounded by an accretion disk. Several obser- 
vations of X-ray binaries have indicated the presence of 
NSs in binary systems in which the NS accretes from 
the less-compact companion (e.g., iBildsten et aT1ll997l ). 
Furthermore, isolated NSs are likely to be surrounded 
by a fallback disk re mnant of the supernova explosion 
(e.g., [Chevalier 198§), while numerical simulations have 
shown that the accretion induced collapse (AIC) of a 
white dwarf may lead to an NS surrounded by a d isk 
with mass up to ~ 0.8M Q (jAbdikamalov et al.ll2uTot) . 

Even if the NS has initially a mass below the maxi- 
mum limit, it can soon become unstable to gravitational 
collapse by accreting matter from the surrounding disk. 
Such AIC to a spinning black hole (BH) could be a pow- 
erful source of gravitational waves (GWs), especially if 
the accreting matter excites non-axisymmetric modes in 
the NS before collapse. To date, no studies of such 
sources have been reported in full general relativity. On 
the other hand, the collapse of NSs (in vacuum) to BHs 
has been the subject of several investigations, including 
the accurate extraction of the G W signals emitted by 
these sources dStark fc Piranlll985HShibata et all 120001: 
Shibatall200l iBaiotti et all I2005allb1: IBaiotti fc Rezzollal 



20061: IBaiotti et all 120071: iGiacomazzo et all 1201 lbO . All 
these simulations have shown that the collapse is essen- 
tially axisymmetric and that no torus is left behind af- 
ter the creation of a spinning BH. Because of this, the 
GW signals emitted by these sources will be rather weak 
and detectable only by advanced LIGO and Virgo if the 
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sources were located in our Galaxy. Third generation 
detectors, such as the planned Einstein Telescope, could 
instead d etect some of these source s up to a distance of 
~ 1 Mpc (|Giacomazzo et alll2TJTlB . 

The AIC of an NS to a BH has also been sug- 
gested to be behind the cen tral engine of short gam ma- 
ray bursts (GRBs e.g. , see iMacFadven et al.l 120051 and 
iDermer k Atovanll2006D . The idea is that an NS in a low- 
mass binary system could accrete enough mass from the 
companion star to produce a spinning BH surrounded by 
an accretion disk massive enough to power a relativistic 
jet and hence a short GRB. 

Here we present the first numerical simulations describ- 
ing the AIC of an NS to BH via fully three-dimensional 
general relativistic simulations. Wc follow the collapse 
of three initial NS configurations, each with and with- 
out an accreting torus, and study the impact of the ac- 
creting material on the dynamics of the collapse and its 
imprint on the GW signal, as well as the post-collapse 
evolution of the disk. Given that this is the first of such 
studies, we consider non-magnetized NSs. We also note 
that, while we do not expect magnetic fields to have a 
strong impact during the collapse (because, even when 
considering highly magnetized NSs, the ratio of magnetic 
to gas pressure is much lower than one), they can have 
an important role after BH formation since they may 
extr act energy from the system and power relativistic 
jets (iRezzolla et alll2011l) . 

Our Letter is organized as follows: in Section [2] wc 
describe the initial NS/torus configurations, and the nu- 
merical setup used to study them. The dynamics of the 
collapse is reported in Section [3J while the GW signal 
emitted during the AIC event is described in Section [U 
followed by a discussion (Section [5]) of possible electro- 
magnetic (EM) counterparts. We finally summarize and 
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Table 1 

Properties of the Initial Models 



Model J c r v /r e W b M R~ e J/W 2 C k, 5 R~ R~^~ t M torus 

(10 15 gem- 3 ) (M e ) (M ) (km) (km) (km) (Mq) 

la 1.91 0.85 1.90 1.73 12.2 0.36 0.21 

lb 1.91 0.85 1.90 1.73 12.2 0.36 0.21 3.63 -0.60 13.6 31.6 0.141 

2a 1.85 0.75 1.98 1.80 13.2 0.47 0.20 

2b 1.85 0.75 1.98 1.80 13.2 0.47 0.20 3.61 -0.74 14.7 32.5 0.143 

3a 1.85 0.65 2.05 1.86 14.4 0.55 0.19 

3b L85 0.65 2.05 1.86 14.4 0.55 0.19 3.60 -0.80 15.5 33.4 0.145 

Note. — The columns refer respectively to the name of the model, the central value p c of the rest-mass density of the NS, the ratio 
r p /r e of the polar to equatorial radius of the NS, the baryonic (M^) and gravitational (M) masses of the NS, the proper equatorial radius 
R e , the ratio J/M 2 of total angular momentum and gravitational mass, and the compactness C = M/R e . The disk is characterized by the 
parameters Iq and S (see the text for description), its inner (Ri n ) and outer (Rout) radius, and its baryonic mass Mtorus- 



conclude in Section [5] 

Throughout this Letter we use a spacelikc signature of 
(—,+,+,+) and a system of units in which c = G = 
Mq = 1 (unless specified otherwise). 

2. INITIAL DATA AND NUMERICAL SETUP 

The initial data are built u sing the cod es RNS 
(iStergioulas k Friedman] 119951: iNozawa et al.l 119981 : 
IStergioulasI I2003D and TORERO (jCorvinol 120091 ? More 
specifically, the uniformly rotating NS models are built 
using the numerical code RNS with a polytropic equation 
of state (EOS), p = Kp r . where p and p arc pressure 
and rest-mass density, respectively, while K = 100 and 
r = 2 are the polytropic constant and the polytropic 
exponent, respectively. All the i nitial models are cho sen 
to be on the unstable branch (Ta kami et al.l 120 111) in 
order to be able to compare the GW signal between a 
collapsing NS with and without a torus. In the former 
case we assume that the torus has already accreted 
enough mass on the NS to make it unstable, and we 
follow its dynamics after that point. The NS collapse 
is triggered by reducing its pressure by 0.1% at the 
beginning of the simulations. 

When we add an accretion disk to the rotating NS, 
we use the code TORERO, which can compute stationary 
solutions of axisymmctric disks ar ound spinning BHs by 
solv ing the equation s described in lDaigne k Font! ()2004l ) 
and iCorvinol ([2009) . In these equations the self-gravity 
of the disk is neglected and hence they are valid only 
when the mass of the disk is sufficiently small compared 
to the mass of the central compact object. We note that, 
even if the torus is built assuming the spacetime of a 
spinning BH, the violation of the Hamiltonian constraint 
introduced by having the torus orbiting around a spin- 
ning NS is not large enough to affect the dynamics of the 
system. This is due to the fact that the tori considered 
here have masses which arc a small fraction of the mass 
of the NS. We have further verified that this violation 
decreases with increasing resolution and that it is negli- 
gible compared to the violations introduced during the 
evolution (i.e., the infinity norm of the Hamiltonian con- 
straint is similar during the collapse between the models 
with and without torus). Furthermore, the initial viola- 
tion produced by adding the torus produces a spurious 
burst of radiation that leaves the grid well before the 
beginning of the collapse and hence it does not contam- 
inate the GW signal computed from these models. The 
tori are built by giving as input parameters to TORERO 
the gravitational mass M and the ratio of angular mo- 



mentum to gravitational mass, J/M 2 , of the NS, and two 
parameters Iq and S which define, respectively, the initial 
value of the constant angular momentum in the disk and 
the depth level reached by the disk inside the potential 
well. A negative value of 5 represents a stable solution, 
while a positive value would give a torus overflowing its 
Roche lobe. 

Table [T] lists the relevant properties of the NSs and tori 
studied here. We have chosen tori with masses ~ 7% of 
the mass of the NS in order to have a non-negligible effect 
on the dynamics. 

The evolution of these NS/torus configurations is 
computed using the fully gener al relativisti c mag- 
netohydrodynamic code Whisky dBaiotti et all [2005at 
IGiacomazzo k Rezzollall2007t iGiacomazzo et alll2011al) . 
which adopts a flux-conservative formulation of the equa- 
tions and high-resolution shock-capturing schemes. All 
the results presented here have been obtained by means 
of the Piecewise Parabolic Method, while the Harten- 
Lax-van Leer-Einfeldt approximate Riemann solver has 
been used to compute the fl uxes (jBaiotti et al.l l2005al : 
IGiacomazzo k Rezz olla 20071 ) . All the simulations were 
performed using an ideal- fluid EOS with T = 2. Mag- 
netic field effects were not included in these simulations 
and they will be the subject of future work. The evo- 
lution of the spacetime was obtained using the Ccatie 
code, a three-dimensional finite-differencing code provid- 
ing the solution of a co nformal traceless for mulation of 
the Einstein equations (jPollney et al.l 120071 ) . Both the 
Ccatie and W hisky codes benefi t of the use of the Ein- 
stein Toolkit (jLoffier et al.ll2012l ). 

The simulations were run using the vertex-centered 
fixed m esh-refinement approa ch provided by the Carpet 
driver (jSchnetter et al.l 120041 ) . We have used seven re- 
finement levels with the finest resolution being 0.1 Mq pa 
0.148 km and the coarsest resolution being 6.8M Q pa 
10.04 km. The finest grid has a radius of lOM© pa 
14.77 km (which is sufficient to completely cover the 
NS), whereas the coarsest grid extends to 998. 4M Q pa 
1474 km. For all the simulations reported here we have 
also used a reflection-symmetry condition across the 
z = Oplane and a 7r-symmetry condition across the x = 
plan<o 

3. DYNAMICS 

3 Stated differently, we evolve only the region {x > 0, z > 0} 
and apply a 180°-rotational-symmetry boundary condition across 
the plane at x = 0. 
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Figure 1. Evolution of the rest-mass density p on the equatorial plane for models 3a (left panels) and 3b (right panels). The top panels 
show the initial condition, the central one the moment of BH formation, and the bottom ones the end of the simulations. In the central 
and bottom panels, a thick black line denotes the position of the apparent horizon of the BH. 
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Figure 2. Evolution of mass (bottom panels) and accretion rate (top panels) for models la and lb (left panel), 2a and 2b (center panel), 
and 3a and 3b (right panel). In all the panels the blue solid line refers to the case without torus (la, 2a, 3a) and the red dashed line to 
the case with the torus (lb, 2b, 3b). 



Figure Q] shows the evolution of the rest-mass density p 
on the equatorial xy plane for models 3a (left panels) and 
3b (right panels). We remind the reader that the only 
difference between the initial data of these two models 
is the fact that we added a torus in the case of model 
3b. The top panels show the initial data, the center ones 
the time of collapse to BH, and the bottom panels the 
end of the simulations. While the dynamics is essen- 
tially axisymmetricQ as it was alread y seen in the case of 
the collapse of NSs with out a disk dBaiotti et al.l [20071 : 
iGiacomazzo et aill2011bl ) , the main difference is that the 
torus resists the gravitational collapse and continues to 
orbit and accrete onto the spinning BH (bottom-right 
panel). Moreover, since the accreting torus triggers oscil- 
lations in the NS density, the collapse is slightly delayed 
in the case in which the NS is surrounded by a disk. This 
happens because the NS oscillations prevent the central 
density from growing exponentially from the very begin- 
ning, since each NS oscillation produces an expansion 
and hence a small decrease of its central density. This 
qualitative dynamics is similar also in the other cases 
we studied (models la/lb and 2a/2b). In the bottom 
panels of Figure [2] we show the evolution of the total 
baryonic mass M&, normalized to its initial value. All 
the cases which included a torus at the beginning (red 
dashed lines) still preserve ~ 7% of the initial mass out- 
side the BH at the end of the simulations. On the other 
hand, all the models without a torus (blue solid lines) 
do not leave any mass outside the BH, as it was already 
seen in previous si mulations of the collapse of an NS to a 
BH in vacuum^! (|Baiotti et al.l l2007t IGiacomazzo et all 
l2011b| ) . Another important quantity is the accretion rate 
Mb after the formation of the BH (top panels of figure[2]). 
We compute this as the integral of the flux of matter that 
crosses the apparent horizon. In all the models that left 
a torus outside the BH, its early-time value is found to 
be - 10~ 2 M Q s" 1 . 

4 In order to assess the impact of the use of 7r-symmetry on our 
results, we also performed some simulations without using it and 
we verified that the dynamics is not affected and that the evolution 
remains axisymmetric. 

5 Note that the presence of magnetic fields alters these con- 
clusions. Simulat ions of magnetized, differentially rotating NSs 
HDuez et al.H2006al lB show that, following the NS collapse, the BH 
is surrounded by a hot torus. 



We note that in our simulations (as well as in previous 
general-rclativistic simulations of NS-NS mergers), the 
accretion is triggered both by numerical viscosity and by 
the spacetime dynamics. In order to assess the contami- 
nation from the numerical component (i.e., the numerical 
viscosity), we performed simulations at different resolu- 
tions, and found that the accretion rate remained at the 
same level. Hence we can safely conclude that the con- 
tribution from numerical dissipation is negligible in our 
results. However note that, in astrophysical disks, turbu- 
lent viscosi ty, likely generated by th e magneto-rotational 
instability (jBalbus fc Hawlevlll991l ). causes angular mo- 
mentum transport and hence accretion. For an a- 
viscosity parameter ^0.1 (iShakura fc Sun vacv 1973), we 
found that the viscous timescale in our torus is to ~ 0.1 s, 
which would yield an even larger post-collapse accretion 
rate M ~ Mt orus /fo ~ 1 M Q s _1 . Since this estimate is 
highly dependent on the disk properties and its viscosity, 
we note that our computed post-collapse Mb should be 
rather considered as a lower limit to what can be achieved 
by the system following an AIC event. 

4. GRAVITATIONAL WAVES 

Figure [3] displays the GW signal emitted by the AIC 
of an NS to a BH in our six models, and in particular 
the dominant I = 2, m = mode computed using th e 
Weyl scalar * 4 (jBaker et al.l 120021 : iLofner et all 120121 ). 
The blue solid lines refer to the cases without a torus, 
while the red dashed lines show the signal emitted when 
the torus is present. Since the collapse is essentially 
axisymmetric, / = 2, m = is the dominant mode, and 
the signal is given by an exponential increase followed 
by the ring-down of the BH. No significant difference 
can be seen in the signal between the models with 
and without a torus, except for the fact that the GWs 
emitted by the models with a torus show some small 
oscillations before the collapse. These are related to the 
NS oscillations triggered by the accretion of the disk, as 
previously mentioned in Section [3] Similarly to the case 
of the collapse of an NS in vacuum, the GW emission 
would be detectable by advanced LIGO/Virgo only if 
the sources were located in our Galaxy. On the other 
hand, with the Einstein Telescope, the signal could be 
detected for sources at distances of up to ~ 1 Mpc. 
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Figure 3. Gravitational wave signal /C)(i = 2,m=0) f° r models la and lb (left panel), 2a and 2b (center panel), and 3a and 3b (right 
panel). All the signals are extracted at r = IOOMq ss 148km. In all the panels the blue solid line refers to the case without torus (la, 2a, 
3a) and the red dashed line to the case with the torus (lb, 2b, 3b). The time is shifted by the time of collapse t C ollapse (defined as the 
moment at which an apparent horizon is formed). 

in an AIC is much shorter than the flare d uration), un- 
less the disk undergoes som e instability (e.g. iPerna et al.l 
[20(31 iProgalT Zhang 2006). 



5. ELECTROMAGNETIC EMISSION 

As described in Section [3l the presence of a rapidly 
accreting torus left around the newly born BH is a dis- 
tinctive feature of the AIC, differentiating it from the 
direct collapse of an NS into a BH, which leaves no mass 
for the BH to accrete from. 

The accretion rates that our AIC events predict (cf. 
Figure [2|), on the order of <~ 10~ 2 Af Q s , are within an 
order of mag nitude of those obtai ned in the merger of 
two NSs fe.g.. iRezzolla et al.ll2010l) . For a typical mass- 
to-energy conversion efficiency ~ 0.1, these AICs have 
luminosities sufficiently high to constitute enticing can- 
didates for short GRBs. In fact, assuming an energy-to- 
7-ray conversion efficiency ~ se yeral tens of perce nt (as 
expected for short GRBs, e.g., IZhang et al.l 120071 ). our 
simulated AIC events would have peak luminosities in 
7-rays on the order of ~ 10 51 erg s _1 , compatible with 
at least a fraction of short GRBs. Clearly, the predicted 
accretion rates are dependent on the torus mass, and a 
more comprehensive exploration will be performed in fu- 
ture work. 

Short GRBs from the AIC of an NS in a binary are 
envisi oned to have the observe r located out of the binary 
plane (jMacFadven et al.ll2005f ). and a jet launched along 
the NS rotati on axis - either as the result of vv anni- 
hilation (e.g., lAlov et ail 120051) or of magneto hydrody- 
namic effects (e.g.. IDrenkhahn fc Spruitl 120021) : further- 
more, they are expected to be more o ften found in early- 
type or globular cluster environments (|Dermer fc Atovanl 
l2006h . A detailed computation of their expected cosmo- 
logical rates, and a comparison with the observed pop- 
ulation of short GRBs, is reserved for future work. In 
particular, we will examine the dependence of the AIC 
accretion rates on the disk mass, and hence assess the 
likelihood of their occurrence in X-ray binaries (where 
disk masses are smaller than the ones cons idered here) 
and f allback disks around isolated NSs (e.g.. IPerna et al.l 
120001 and references therein). Finally, within the con- 
nection short GRBs/AIC events, we further n ote that 
the observed X-ra y flares (iMargutti et al.ll201lf) and ex- 
tended emissions (|Norris fc Bonnelll 120061 ). lasting up to 
several hundreds of seconds, cannot be directly explained 
as the result of activity from the central engine (which 



6. CONCLUSIONS 

We have performed the first three-dimensional, fully 
general relativistic simulation of the AIC of an NS to 
a BH; we have examined the GW signal expected from 
such events, and computed the accretion rates onto the 
newly formed BH. We find that the main component 
(I = 2, to = mode) of the GW signal is very similar 
to the one from the collapse of an NS into a BH in vac- 
uum. However, and most importantly, while the direct 
collapse of a non-magnetized NS onto a BH leaves no 
mass behind (hence making unlikely the presence of any 
EM detectable counterpart), an AIC is followed by rapid 
accretion of the disk onto the newly formed BH. Our 
computed accretion rates, comparable to those obtained 
in simulations of NS-NS mergers, make AIC events po- 
tenti al engines for short G RBs. NSs in X-ray binaries 
(e.g., iBildsten et al.l 119971), is olated NSs accreting from 
fallback disks ( Chevalierlll989D, and NSs produce d by the 
AIC of a white dwarf ( Abdikamalov et al.ll2010l ) are as- 
trophysical scenarios which might host the AIC of an NS 
onto a BH. 
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